Gravitational Cherenkov losses in MOND theories 
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Survival of high-energy cosmic rays (HECRs) against gravitational Cherenkov losses is shown 
not to cast strong constraints on MOND theories that are compatible with general relativity (GR): 
theories that coincide with GR in the high-acceleration limit. The energy-loss rate, E, is shown 
to be many orders smaller than those derived in the literature for theories with no extra scale. 
The gravitational acceleration produced by a HECR in its vicinity is much higher than the MOND 
acceleration ao- So, modification to GR, which underlies E, enters only beyond the MOND radius 
of the particle, within which GR holds sway: r M = (Gp/cao) 1 ^ 2 . The spectral cutoff, which enters E 
quadratically, is thus r^/, not the particle's, much larger, de Broglie wavenumber: k dB = p/h. Thus, 
E is smaller than published rates, which use k dB , by a factor (r M k dB ) 2 « 10 39 (cp/3 x 10 n Gev) 3 



With r M as cutoff, the distance a HECR can travel without major losses is q£ M , where £a 



is the MOND length, and q is a dimensionless function of parameters of the problem. Since l M is 
~ 2-7T times the Hubble distance, survival of HECRs does not strongly constrain GR-compatible, 
MOND theories. Such theories also easily satisfy existing preferred-frame limits, inasmuch as these 
limits are gotten in high-acceleration systems. I exemplify the results with MOND adaptations of 
Einstein- Aether theories. 
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PACS numbers: 04.50.Kd, 04.30.- 



I. INTRODUCTION 



98.70.Sa 



Several modifications of general relativity (GR) have 
been considered in recent years that involve nonluminal 
gravitational wave propagation. Depending on the the- 
ory, and on the choice of its parameters, these nonluminal 
waves can propagate faster or slower than light. Exam- 
ples of such theories are Einstein- Aether (EA) theories 
(e.g., [1], and references therein). Superluminal propaga- 
tion is usually frowned upon, because it connotes closed 
timelike curves with all their problematics, but dissent- 
ing opinions have also been heard (e.g., [2-5]). Sublu- 
minality, on the other hand, may be in conflict with the 
observation of high- energy cosmic rays (HECRs) that are 
thought to arrive to earth from astrophysical distances. 
Caves [6] , and, more recently, Moore and Nelson [7] , have 
pointed out that unless the subluminal speed is very near 
the speed of light, such HECRs, which themselves move 
very near the speed of light, lose their energy very ef- 
ficiently by gravitational Cherenkov radiation (GCR) of 
the subluminal waves. Elliott et al. [8] then discussed 
the constraints that such considerations cast on the pa- 
rameters of EA theories. 

The relativistic formulations of MOND proposed to 
date, in the form of modified-gravity, involve gravi- 
tational degrees of freedom in addition to the met- 
ric. 1 These theories include TeVeS and its predecessors 
([9, 10], reviewed in [11]), adaptations of EA theories 
to MOND [12], and BIMOND [13]. Wave propagation 
in BIMOND has not been studied yet, but in the other 
MOND theories there are nonluminal modes. So, either 
parameters can be chosen so that these waves are always 



superluminal, or else, these parameter need to be sub- 
jected to the constraints from HECR survival. 

I point out here that results as in refs. [7, 8] cannot be 
applied to constrain MOND theories, as has been stated 
before, without major corrections that greatly weaken 
the resulting constraints. 

In section II, I show the need for a new derivation of 
the loss rate in MOND theories, and obtain the new esti- 
mate. In section III, I exemplify the results with MOND 
adaptations of EA theories. 



II. MOND ADJUSTED LOSS RATES 

Detailed calculations of the GCR energy-loss rates by 
HECRs have been described in [6-8] , for modifications of 
GR with no extra scale. Cutting through the details, we 
may note that as the electromagnetic Cherenkov emis- 
sion is proportional to the square of the particle charge, 
so the GCR rate is proportional to G{^M) 2 = G(p/c) 2 , 
where G is Newton's constant, and p and 7 are the par- 
ticle's momentum and Lorentz factor. Then, on dimen- 
sional grounds, the energy radiated per unit time per unit 
wavevector, fc, is 



d 2 E _ SG P 2 k 
dkdt c 



(1) 



where S is some dimensionless function of parameters 
of the theory, such as the wave speed for the particular 
mode (in units of c), of hk/p, and of the particle's j3 = 
v/c. To get the total energy-loss rate, this expression 
is integrated over k, up to some cutoff. In [7, 8] this is 
k dB = p/h. The resulting expression for the loss distance 
is 



1 This is not necessarily the case in "modified inertia" versions. 



(2) 
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where Q is a dimensionless quantity that depends on di- 
mcnsionless parameters of the theory, the wave mode un- 
der consideration and on the speeds of the particle and 
the wave in units of c. The coefficient Q depends also on 
the internal structure of the particle (e.g., partonic struc- 
ture of a hadron) . This enters because the radiation pro- 
cess assumes coherence over the radiating particle, but 
for a high-energy hadron, the length k~^ B already probes 
its internal structure; so coherence applies only to the 
smaller constituents. In [7] this is taken into account by 
considering radiation from partons that each carry 0.1 of 
the hadron momentum. This increases Di oss by a factor 
10 3 over that for the whole hadron; in [8] the structure 
is treated in more detail. Equipped with expression (2), 
one uses lower limits on the travel distance, D, of high 
p cosmic rays, to impose Di oss > D, thus placing upper 
limits on Q, and hence on the parameters of the theory. 

Here I note that the situation in relativistic MOND 
modifications of GR is very different, because they intro- 
duce a new dimensioned constant: the MOND acceler- 
ation a w 1.2 x 10~ 8 cms~ 2 . I consider only MOND 
theories that coincide with GR in the limit clq — > 
called here GR-compatible, MOND theories. This can 
be straightforwardly implemented in the MOND adapta- 
tions of EA theories, in BIMOND, and in the theories 
discussed in [14]. It is not quite the case for TeVeS, 
where a departure from GR remains even in the limit 
(see, e.g., [9, 15]). In GR-compatible theories we expect 
an attenuation of the GCR in circumstances where the 
gravitational accelerations are larger than do- For exam- 
ple, a particle traversing the inner solar system, where the 
gravitational acceleration due to the sun is much higher 
than a , would emit strongly attenuated GCR, except, 
perhaps, at wavelengths longer than the MOND radius 
of the sun, r„, = (M G/a o ) 1/2 « 10 17 cm. 

The shorter-wavelength "anomalous" waves, whether 
subluminal or superluminal, have a much suppressed in- 
teraction with matter in high-acceleration backgrounds, 
such as the inner solar system, or the vicinity of close 
binary pulsars. So, their detectibility, and emissivity, in 
such systems is greatly suppressed. 

The high-acceleration regime is also very pertinent for 
hadrons: The gravitational acceleration produced by a 
hadron of mass M, at rest, at a distance h/Mc, is w 
2.6 x 10 4 (M/M p ) 3 a , where M p is the proton mass. For 
a high-energy hadron of Lorentz factor 7, the near-field, 
gravitational acceleration it produces at the distance k dB 
is a factor -f 3 larger, which is many orders of magnitude 
larger than ao- 2 

It is impracticable to calculate exactly the GCR loss 
rate in any particular MOND theory, let alone give a gen- 
eral result. We can, however, make an estimate based 
on the following argumentation: I surmise that for GR- 



2 I am speaking here of the near field of the hadron, not the radi- 
ation field. 



compatible, MOND theories, there are no Cherenkov ef- 
fects within the MOND radius of the particle, where GR 
holds sway: 

r M = (Gp/ca Q ) 1/2 « 3 x l(T 12 (cp/Gcv) 1/2 cm. (3) 

Effectively, the particle carries around it a bubble of ra- 
dius r M , devoid of the subluminal waves. Waves with 
wavenumbers larger than r^ 1 arc then not produced co- 
herently. I also assume that in the MOND regime, far 
beyond r M , the theory has no further scale, and so we 
can use the results from the literature for the gravita- 
tional Cherenkov spectrum for such theories. This means 
that we can adopt the Cherenkov spectrum discussed 
above, but cut it off above the wavenumber ~ , in- 
stead of at kdB- The MOND radius is many orders of 
magnitude larger than kj B : r M k dB w 160(cp/Gev) 3 / 2 . 
In the MOND context the rates appearing in the lit- 
erature should thus be reduced, nominally, by a factor 
~ (r M k dB ) 2 w 10 39 (cp/3 x 10 n Gev) 3 , normalized to the 
momentum used in the constraints in [7, 8] . 

More precisely, to account for MOND effects we may 
approximately include in expression (1), a factor 1 — 
fi(a/ao), where, (for a > a ) a/a — {kr M ) 2 measures the 
acceleration at the radius k~ l probed by wavenumber k, 
and /j, is a MOND interpolating function appropriate for 
the circumstances. So the factor 1 — fi measures the de- 
parture from GR at a distance The integral over 
k now converges if 1 — /j,(x) vanishes at high x, faster 
than x^ 1 , which is known to be required, anyhow, from 
solar-system constraints. The MOND factor effectively 
cuts off the integral at k = r^ 1 . So, integrating the 
Cherenkov radiation spectrum, Eq.(l), up to the MOND 
cutoff, instead of the de Broglie wavenumber, expression 
(2) is replaced by 3 

D loss = qt M . (4) 

Here, £ M = c 2 /a is the MOND length, and q, which 
has a similar role to that of Q^ 1 in Eq.(2), depends on 
parameters of the theory and on the wave and particle 
velocities (in units of c). 4 Observations that place a lower 
limit on Di oss then place lower limits on q. As has been 
noted before in other contexts, £ M « 2itDh, where Dh is 
the Hubble distance; so the lower limits on q are of order 
unity, even if wc require survival over Dh- 

For low-rest-mass particles, such as electrons, photons, 
and neutrinos, whose de Broglie wavelength can be larger 
than their MOND radius, high-acceleration effects enter 
only if r M k dB > 1, which is the case for energies E > 
E M = (c 4 ft 2 o /G) 1/3 = E p (t p a / 'c) 1/3 - 35Mev, where 
E p and t p are the Planck energy and time, respectively. 
Then our expression above for Di oss applies. For lower 



3 Since the fc-cutoff is now much smaller than p/H, we ignore the 
dependence of S on kh/p, which can be set to in the argument. 

4 The factor q insures, for example, that Di oss is infinite if the 
particle speed is below the wave speed. 
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energies, r M k dB < 1, and it is k dB that sets the cutoff. 
Then, the expression for D loss is multiplied by (E M /E) 3 . 

Since the MOND radius is much larger than the size 
of a hadron, the particle radiates coherently, and there 
is no need to account for its internal (parton) structure. 
So, after correcting the rates for scaleless theories by a 
factor (r M k dB ) 2 , to make them fit for MOND theories, 
we also have to correct back the partonic effects on the 
rates. For example, in [7], Q defined in Eq.(2) is Q ~ 
10~ 3 (n — l) 2 , where n is the ratio of c to the wave speed 
(the particle speed is assumed much nearer to c) , and the 
factor 10 -3 assumes that the relevant partons carry 0.1 
of the hadron momentum p. In the MOND context we 
would have, instead, q^ 1 ~ (n — l) 2 , where the factor 
10~ 3 has been removed. So even for travel distances of 
the Hubble length, the bound on n — 1 is only of order 
1. (The expression for Q in [7] assumes n — 1 <C 1; so a 
more careful treatment is needed in general.) 

The effect I propose here is analogous to that calcu- 
lated for the electromagnetic Cherenkov case, when a 
charged particle moves along the axis of a vacuum cylin- 
der of radius R (with refractive index n = 1) surrounded 
by a medium with n > 1. Since the medium with n > 1 
is kept a distance R from the particle, the wavenumber 
spectrum is, indeed, found to be sharply cutoff above 
~ R^ 1 (see e.g., [16]). This happens also for a particle 
moving in vacuum, parallel to the interface of a half- 
space dielectric, with n > 1, a distance R away. Per- 
haps a better analogy is when the particle, which moves 
through a medium with n > 1, induces around it, and 
carries with it, a region of radius R (in the medium's rest 
frame) where n = 1 to a high accuracy. My result is also 
akin to the effect of a finite radius, R, of the emitting 
charge (e.g., [17, 18]), where the radiation spectrum is 
also sharply cut off for kR > 1. 

There is, clearly, a more general lesson to be learnt 
from our analysis. For example, similar corrections would 
apply to other GR-compatible modifications involving a 
new scale, such as modifications of GR at large distances. 

Note, in passing, that constraints from solar-system 
and binary-pulsar limits on preferred-frame parameters, 
are easily satisfied by GR-compatible, MOND theories 
(not by TeVeS, [15]). These limits, inasmuch as they 
pertain to experiments in high-acceleration regions, test 
such theories where they differ from GR by only ~ 1 — 
(i(a/ao), which can be extremely small. Note also, that 
in the MOND limit of such theories (a — >• oo, with Gao 
fixed), the post-Newtonian expansion is not applicable, 
since the nonrelativistic (NR) limit of the theory is then 
not Newtonian, contrary to the basic premise of the post- 
Newtonian expansion. 



III. EXAMPLE: MOND ADAPTATION OF 
EINSTEIN-AETHER THEORIES 

Consider, as a specific example, the MOND formula- 
tions that were adapted from EA theories [12]. In EA 



theories ([1], and references therein), gravity is carried 
by a metric, g ilVl as well as a vector field, A a . One adds 
to the standard Einstein-Hilbert Lagrangian for the met- 
ric, the terms (in this section I take c = 1) 

C ( A >9) = ^Lk + £l, (5) 



JC = a - 2 /C^ ;Q A%, (6) 

= c ig af) g ia + c 2 5*5P + c 3 5^ + c 4 A a APg ja , (7) 

and £i is a Lagrange multiplier term that forces the 
vector to be of unit length. 

Elliott et al. [8] considered three constraints on the 
coefficients, resulting from GCR losses to the different 
modes. 

These EA theories have been adapted to give MOND 
phenomenology in the NR regime [12], by modifying the 
first term in the added Lagrangian to 

C M {A,g) = ^F{lC). (8) 

Zlosnik et al. [12] limited themselves to theories with 
c 4 = 0. They showed that to get the correct NR MOND 
phenomenology we have to have T{K) « K. + r\K?l' 1 for 
K, <C 1, and c\ — —2 (with their sign convention for the 
Ci, which is opposite to that of [8] and of others). 5,6 In 
addition, if we want the NR limit of the theory to tend 
to Newtonian gravity for ao — ¥ 0, as it should, we need 
to have T'(K —¥ oo) — »■ 0. In fact, F'{K) plays the role 
of 1 — /i in the discussion leading to Eq.(4). To make 
the theory GR compatible, I require, to be on the safe 
side, that T is bounded by a constant £. Matter (and the 
metric) then decouples from the vector field in the limit, 
and is governed by GR. 7 (If T — > £ ^ 0, we get GR with 
a cosmological constant ~ C a o-) Here, as in [12], G is 
Newton's constant. If we replace it with G' ^ G, the NR 
limit can still be made a good MOND theory. To then 
get the Newtonian limit for a — > 0, we have to have 
T'{K, — > oo) — > S 0, with an appropriate S. But then, 
the relativistic theory seems not to be GR compatible. 
A similar situation occurs in BIMOND [13]. 8 



5 This makes the speed of one of the modes-the trace aether-metric 
mode-vanish, assuming no background acceleration field, W<j>b- 
If our system is in such a field, e.g., in the vicinity of galactic 
systems, (u/c) 2 of the different waves may suffer additive correc- 
tions of order |V<fe|/ao, when |V0(,|/arj < 1. 

6 A constant factor a in front of K on the right-hand side, as 
appears below Eq.(10) of [12], is absorbed in the Cj, and is taken 
a = l. 

7 Thus, vector waves will be very difficult to detect in the inner 
solar system, if they have A; -1 < r M ,o- 

8 Interestingly, MOND adaptations of EA are constrained versions 
of the BIMOND class: If we force the two metrics g^" and g^" , 
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It is impracticable to get an exact solution of the full 
problem of GCR, especially as the problem involves the 
full acceleration range, as we saw above. As a very rea- 
sonable approximation I proceed as follows. In the field 
equations, the Cj, when they are outside the argument 
of J 7 , appear always multiplied by J"'(/C). So, in high- 
acceleration environment, where J"'(/C) — > 0, we are in 
a situation that can be reproduced by taking the limit 
Cj — > 0, with their ratios fixed, in a standard (scale free) 
EA theory. In this limit, all the loss rates calculated 
in [8] vanish, as expected. This implies that the GCR 
losses are, indeed, cut off for waves with k > k M — v M . 
(As explained earlier, the cutoff is not sharp, but since 
the Cherenkov spectrum increases only as fc, the cutoff is 
effectively at k M .) 

Next, we note that in the MOND regime-where K. <C 1, 
and so J-(IC) « /C-we get an EA theory, as studied by 
[8]. 

I thus approximate the rates in this MOND adapta- 
tion of EA theories, by making two corrections to the 
expressions for the rates in [8]: First I replace their cut- 
off wavenumber by r^ 1 , as described in section II. Then I 
increase the rate by a factor of 10 3 to account for the fact 
that with the smaller cutoff, the hadron radiates coher- 
ently, and not as the sum of its partonic contributions. 
This factor corrects back the partonic correction of [7]. 9 

The constraints derived in [8] are in the form of upper 
limits on three functions of the Cj, fk{ci) 1 k = 1, 2, 3, with 
coefficients of order 1 (one for each mode), assuming that 
HECRs travel a distance D, of at least lOkpc without 

which in BIMOND arc independent degrees of freedom, to be 
related by = g"" + sA^A v (s < 1), with, furthermore, 

g^A^A" = -1, we get the EA-MOND class of theories, with 
g^v and A^ as degrees of freedom. For example, with these con- 
straints, the BIMOND scalar arguments T, T, shown in [13] 
to lead to particularly simple theories, give EA-MOND versions 
with ci = C4 = — C3, while C2 is free and introduced by scalars 
of the form g^^C^C" , and its image under metric interchange, 



major losses. These functions vanish when Ci = 0, and 
are, generically, of order 1 when the Ci are of order 1. 
Only the forms of ft for a <C 1 are given, because the 
constraints in [8] require this; but this is not the case 
for MOND theories. The constraints in [8] for standard 
EA theories [their eqs.(3.13) (squared here), (4.7), and 
(4.15)] are: 10 

fk(ci) < M-D/lOMpc)- 1 , (9) 

where b l = 2.5 x KT 34 , b 2 = 7 x 10" 35 , b 3 = ICT 33 . 
With the MOND modification we get instead &i = 250, 
62 = 70, 63 = 10 3 . So, even if the coefficients are of 
order I, HECRs can propagate a Hubble distance before 
they lose much of their energy, in line with the result in 
section II. 

In addition to their constraints from GCR, [8] consider 
the emission of a double scalar mediated by a graviton. 
This process will not concern us here, as it involves some 
more delicate processes and assumptions. In any event, 
even if this constraint remains, it is much less difficult to 
satisfy, as it is only one constraint on the ratio of the q , 
not on their absolute values, as for the Cherenkov con- 
straints. It is probably also subject to major relaxation 
in the MOND context, via the mechanism I discuss here. 
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terms which retain the simplicity. 
9 Elliott et al. [8] say that their more detailed treatment of the 
structure gives a somewhat stronger constraint than that of [7]; 
this means that our correcting back in this way is conservative. 
10 Here I normalize the constraints of [8] to a travel distance of 
D = lOMpc, instead of their less constraining lOkpc. 
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